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Phase  Correction  for  Sun-Reflecting  Spherical  Satellite 

Erwin  Schmid 

Geodetic  Research  and  Development  Laboratory 

National  Oceanic  Survey 

abstract.  Correction  formulas  are  developed  that  convert  the  ground-based  camera  measure- 
ments of  the  direction  to  the  center  of  the  light  source  on  a  balloon-type,  spherical  satellite 
to  the  corresponding  direction  to  the  geometrical  center  of  the  satellite.  The  correction  is 
necessary  because,  in  the  case  of  a  diffusively  reflecting  satellite,  the  photographed  satellite 
image  refers  to  the  visible  sun-illuminated  portion  of  the  satellite  surface  and,  for  specular 
reflection,  to  the  location  of  the  sun's  image  on  the  balloon.  The  correction  is  small  but,  as  a 
computable  bias,  is  incorporated  in  the  mathematical  model  of  the  Geometric  Satellite  Trian- 
gulation  World  Net  Program. 


INTRODUCTION 

In  the  Geometric  Satellite  Triangulation  World  Net 
Program  (for  which  field  operations  were  concluded 
late  in  1970)  and  in  the  North  American  Densification 
Net  Program  (now  in  progress),  sunlight-reflective 
balloon  satellites,  100  ft  in  diameter,  such  as  the 
NASA-launched  Echo  I,  Echo  II,  and,  at  present, 
PAGEOS  (Passive  Geodetic  Satellite),  are  photo- 
graphed simultaneously  from  two  or  more  ground  sites 
against  the  background  of  surrounding  stars  for  the 
purpose  of  locating  these  camera  stations  in  the  three- 
dimensional  static  flat  space  defined  in  part  by  the 
right  ascension-declination  system  of  metric  astronomy. 
An  exposition  of  the  method  in  detail  is  being  pre- 
pared for  publication  by  the  Office  of  the  National 
Geodetic  Survey,  a  component  of  NOAA's  National 
Ocean  Survey. 

The  light  energy  of  the  sun  reflected  from  the  mov- 
ing satellite  creates  a  track  on  the  camera  plate  which 
is  chopped  by  camera  shutters  into  a  series  of  point- 
like images,  each  of  which  represents  a  position  of  the 
satellite  in  space  at  an  instant  determined  by  the  as- 
sociated electronic  timing  system.  The  measured  plate 
coordinates  of  the  centroid  of  such  an  image  are  trans- 
formed into  space  coordinates  of  the  light  source  (at 
the  satellite)  for  the  specific  image.  It  is  the  purpose 
of  this  report  to  find  the  center  of  this  light  source  on 
the  satellite  and  hence  to  compute  a  correction  that 
displaces  it  to  the  geometrical  center  of  the  balloon  as 
a  common  target  for  all  stations.  If,  as  in  the  case  of 


Echo  II,  the  portion  of  the  balloon  surface  illuminated 
by  the  sun  and  facing  the  camera  site  reflects  the  sun- 
light diffusely,  the  light  energy  sensitizing  the  plate 
comes  more  or  less  uniformly  from  all  parts  of  what 
would  appear  to  the  eye,  if  sufficiently  magnified,  as  a 
crescent  analogous  to  the  moon's  appearance.  The 
specular  reflective  surfaces  of  Echo  I  and  PAGEOS 
concentrate  the  source  of  light  energy  in  a  point  on 
the  balloon  surface.  Pertinent  corrections  are  devel- 
oped separately. 

These  corrections  are  small,  but  they  represent  a 
computable  bias  and  are  incorporated  in  the  mathe- 
matical model  of  the  Geometric  Satellite  Triangulation 
World  Net  Program. 

DIFFUSIVE  REFLECTION 

The  portion  of  a  satellite's  sun-illuminated  hemi- 
sphere facing  an  observer  changes  continuously  with 
time  in  the  same  manner  and  for  the  same  reason  that 
the  moon  goes  through  its  phases.  Figure  1  shows  the 
spherical  satellite  with  radius  p  and  center  P.  If  the 
direction  to  the  sun  points  to  S  on  the  sphere,  then  the 
illuminated  half  is  bounded  by  the  great  circle  DEC 
whose  pole  is  S.  The  direction  to  the  observer  is  E, 
and,  at  great  distances,  his  outline  of  the  satellite  is 
the  great  circle  CAD  for  which  E  is  the  pole.  This 
aspect  from  the  earth  is  shown  in  figure  2,  which  is 
figure  1  viewed  from  below.  The  visible  illuminated 
portion  of  the  satellite  consists  of  the  spherical  triangle 


Figure  1. — Phase  angle 
of  spherical  satellite. 


DAB  in  figure  1  and  its  symmetrical  counterpart  tri- 
angle CAB. 

The  angle  y  at  the  satellite,  between  the  direction  to 
the  sun  and  to  the  observer,  is  defined  as  the  phase 
angle.  Its  range  is  from  0°  to  180°,  producing  at  these 
limits  full  moon  and  new  moon,  respectively.  The  phase 
itself  is  defined  with  respect  to  the  projection  of  the 
illuminated  portion  onto  the  plane  ACD  as  depicted  in 
figure  2;  it  is  directly  proportional  to  the  amount  of 
illumination  reaching  the  earth  (its  relative  magni- 
tude) and  is  equal  to  the  ratio  AH/2AP  of  these  two 
segments  in  figures  1  and  2. 

From  figure  1,  we  have 

PH  =  PB  cos  Z  BPH  =  p  cos  y. 
Hence,  from  figure  2, 


Phase  = 


AH        AP  +  PH 


2AP 


2AP 


p  -+-  p  cos  7       1  +  cos  7 


2p  2 

Assuming  the  sun's  light  to  be  diffused  in  all  direc- 
tions from  the  satellite's  surface,  the  source  of  the  light 
energy  creating  the  satellite  image  may  be  postulated 


The  satellite  image 
centers  on  P1  which 
I  ies  on  vector    PA 
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Figure  2. — Appearance  from 
earth  of  "moon"  phase. 


to  originate  from  a  point  P'  at  the  midpoint  of  segment 
AH.  Because  the  phase  is  generally  different  for  the 
several  camera  stations,  in  addition  to  being  variable 
with  time,  it  is  necessary  to  transform  the  point  P'  to 
a  point  on  the  satellite  invariant  with  respect  to  phase, 
such  as  its  center  P,  to  achieve  coincidence  in  space 
and  time  for  each  of  the  images  photographed  from 
different  locations.  The  purpose  here  is  to  evaluate  the 
effect  of  the  shift  of  the  observation  from  its  source 
P'  to  the  geometrical  center  P  of  the  satellite,  in  other 
words,  the  effect  of  the  displacement  vector  P'P  on  the 
position  vector  camera-satellite.  It  should  be  noted 
that  this  increment  is  sufficiently  small  to  be  treated  as 
a  first-order  differential;  except  for  sign,  it  is  imma- 
terial whether  this  correction  is  added  to  the  direction 
to  P'  or  to  the  direction  to  P. 

The  direction  cosines  of  a  line  in  space  relative  to 
a  specified  Cartesian  coordinate  system  (XYZ  or 
XiX2X3)  are  the  cosines  of  the  three  angles  which  the 
line  makes  with  the  three  axes  in  the  specified  order. 
The  direction  cosines  are  also  the  three  projections 
(inner  or  dot  product)  onto  these  axes  of  a  unit  vector 
associated  with  this  line,  that  is,  the  Cartesian  com- 
ponents of  the  vector.  Figure  3  shows  still  a  third  ap- 
proach which  is  convenient  for  converting  astronom- 
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i-e.  -  1 

*1  =  cos  b  sin  oC. 

l>2-  cos  o  cos^C 

■^3  =  sin  & 


Figure  3. — Direction  cosines  as  Cartesian  coordinates  on  unit 
sphere. 

ical  (spherical)  coordinates  to  vector  components  or 
Cartesian  coordinates.  If  a  sphere  of  radius  one  is 
assumed  to  be  drawn  around  the  origin  of  the  Carte- 
sian coordinate  system,  then  any  radius  vector  of  the 
sphere  will  be  a  unit  vector,  and  conversely.  Thus,  the 
XYZ  coordinates  of  the  end  point  of  a  radius  vector, 
such  as  1,  will  be  the  direction  cosines  l\,  /o,  h  of  the 
corresponding  line  in  space  as  well  as  the  components 
of  the  vector. 

One  small  difference  in  these  different  approaches 
may  be  noted:  In  ordinary  analytical  geometry,  the 
ambiguity  arising  from  the  three  direction  angles  and 
their  supplements  is  left  unresolved  so  that  a  set  of 
three  direction  cosines  is  interchangeable  with  the  set 
having  the  opposite  sign.  This  is  not  permissible  in  the 
vectorial  approach. 

The  Cartesian  coordinate  system  of  figure  3  is  a 
left-hand,  geocentric,  inertial  system;  that  is,  the 
origin  is  at  the  center  of  the  earth,  the  Z-axis  points 
north,  and  the  Y-axis  is  directed  toward  the  vernal 
equinox.  In  this  system,  the  apparent  right  ascension 
and  declination  of  an  object  in  space  are,  except  for 
parallax,  the  angles  a  and  8  indicated  in  the  figure. 

For  an  object  as  distant  as  the  sun,  a  and  S  as  seen 
from  the  satellite  are  sufficiently  close  for  the  purpose 
(maximum  parallax ~  10  seconds  of  arc)  to  use  the 
geocentric  coordinates  a©,  S©  directly  from  the  ephe- 
meris.  Using  these  values,  it  follows  from  figure  3  that 
the  direction  cosines  of  the  line  earth-sun  and,  specifi- 
cally, the  components  of  the  unit  vector  in  the  direc- 
tion from  earth  to  sun  in  the  geocentric  system  or  in 
any  near-earth,  space-parallel  system,  are 


,/i  =  cos  5©  sin  a© 
h  —  cos  5©  cos  a© 

J3  =  sin  5©.  (1) 

The  direction  from  the  camera  station  to  the  satellite 
is  normally  computed  in  terms  of  azimuth  and  eleva- 
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tion  with  respect  to  a  local  horizonal  system.  The 
geodetic  coordinates  of  the  station  are  used  to  convert 
the  azimuth  and  elevation  into  corresponding  right 
ascension  and  declination.  The  latter  are  coordinates 
in  the  spherical  reference  system  that  corresponds  di- 
rectly to  the  geocentric  Cartesian  system. 

Thus,  if  this  apparent  right  ascension  and  declina- 
tion are  a  and  8,  respectively,  then  the  unit  vector  m 
in  the  direction  from  the  camera  to  the  satellite  (direc- 
tion EP  in  fig.  1)  has  components 


,m,i  =  cos  8  sin  a 

,  m2  =  cos  8  cos  a 

Km3  =  sin  5;  (2) 

and  the  unit  vector  in  the  opposite  direction,  that  is, 
from  P  to  E,  has  the  same  numerical  components  with 
the  opposite  sign.  In  other  words,  the  unit  vector  in  the 
direction  of  the  satellite-camera  is  -m. 

The  cosine  of  the  phase  angle  y  is  therefore  -nvl  or 

cos  7  =  —  {hrrix  +  Z2m2  +  Z3m3).  (3) 

Let  the  unit  vector  in  the  direction  PP'  of  figure  1 
be  designated  n.  Because  n  lies  in  the  plane  of  the  unit 
vector  to  the  sun  1  and  the  unit  vector  to  the  observer 
— m,  it  is  a  linear  combination  of  these,  or 

n  =  Al  -f-  fim 

where  A,  /a  are  undetermined  scalars. 

The  scalar  product  of  this  last  equation  with  in  and 
1  gives 

rn'n  =  0  =  Al'm  +  /*m*m  =  —A  cos  y  -f-  ju. 

l'n  =  cos  (  -  —  y  )  =  sin  7 

=  AM  +  /xm*l  =  A— ix  cos  y. 

From  these  two  equations, 

!0  =  —  X  cos  7  +  m 
sin  7  =  X  — ju  cos  7, 
we  find  A  =  esc  y,  ft,  =  cot  y  so  that 

n  =  (esc  y)l  +  (cot  y)m.  (a\ 

This  result  is  also  readily  apparent  from  figure  1  as 

PA  =  PS  +  S'A 

where  S'  is  the  intersection  of  the  line  PS  extended 
with  the  tangent  to  the  circle  at  A.  Scaling  the  figure 
down  from  radius  p  to  radius  unity,  the  directed  dis- 
tances PA,  PS,  PE  become  the  unit  vectors  n,  1,  —  m, 
respectively.   From  triangle  PAS',  the  length   of  PS* 


is  esc  y  and  the  length  of  S'A  is  cot  y.  Therefore,  we 
have 

PS'  =  esc  yl 
and  S'A  =  cot  ym 

because  the  direction  of  S'A  is  opposite  to  PE.  The 
above  vector  equation  is  therefore  the  equation  (4) . 
The  components  of  n  from  equation  (4)  are 


li  +  mi  cos  y 

Hi 

sin  y 

h  +  m2  cos  y 

7l2 

sin  y 

Z3  +  m3  cos  y 

sin  y 


(5) 


The  length  of  the  segment  PP'  is,  from  figure  1  or  2, 


\PP'\  =  AP  -=:AH 


=  p  — g(p  +  p cos  y) 

=  |(1  -cosy), 


(6) 


so  that  the  vector  PP'  =  -(1  —cos  y)  n 

and  vector  P'P  =  ^(cos y  —  l)n.  (7) 

The  unit  vector  m  in  the  direction  of  the  satellite  is 
given  with  equations  (2)  as  a  function  of  right  ascen- 
sion and  declination.  We  consider  now  the  effect  on 
this  a  and  8  of  a  differential  displacement  dm  of  the 
vector.  From  equations  (2),  it  follows  that 

idmi  =  cos  a  cos  5  da  —  sin  a  sin  b  db 
\  dmz  =  —  sin  a  cos  5  da  —  cos  a  sin  5  db 

\dm3  =  cos  b  db.  (8) 

Because  the  components  of  a  unit  vector  are  func- 
tionally related,  only  two  of  these  equations  are  needed 
to  solve  for  da  and  d8.  We  choose  the  first  two  for  this 
purpose  to  obtain  a  more  symmetrical  result  and  get 


7  cos  a  7 

I  eta  =  rami 

I  cos  b 


sin  a 
cos  b 


dm* 


,.  sin  a  7  cos  a  7 

db  =  -    —. — rami  —  — : — rrfm2. 


sin  b 


sin  b 


(9) 


From  the  relation  rax2  +  m22  +  m32  =  1,  it  follows 
that  midmi  +  m2dm2  +  mzdmz  =  0. 

The  left  side  of  this  equation  can  be  interpreted  as 
the  inner  product  of  the  vectors  m  and  dm.  Hence,  the 


differential  increment  of  a  unit  vector  or,  for  that 
matter,  of  any  vector  of  constant  length  is  in  a  plane 
normal  to  the  vector. 

The  phase  correction  vector  P'P  satisfies  this  condi- 
tion and  is  sufficiently  small  to  be  treated  as  a  differen- 
tial. If  the  distance  to  the  satellite  is  the  scalar  D,  then 
the  position  vector  is  Dm  and  its  differential  is  Ddm. 
Setting  this  differential  equal  to  the  vector  P'P  from 
equation  (7),  we  obtain 

dm  =M)  (cos?  _  !)n- 


From  equations  (5),  the  components  of  this  vector 


are 


p(cosy  — 1)  „  . 

=      2D  sin  y      &  + «*«»*) 

p(cosy  —  1)  ,,  x 

dm2  =      nT^    .        -  (l2  +  m2  cos  y) 


2D  sin  y 


p(cos  y  —1)  ,-,     ,  .         .„ .. 

'      *  =      2D  sin  y      (*3  +  m%  C0S  7)'       (10) 

which,  when  substituted  in  equations  (9),  give 

,  p(cosy  —  1)    r  . 

da  =  2D  sin  y  cos  5  [c0S  "{h  +  ^  C0S  7) 

—  sin  a(U  +  m2  cos  y)] 
,.        p(cos  y  —  1)    r  .       ,,  . 

d5  =  2Dsinysin5[sma(Zl  +  miCOS7) 

+  cos  a(h  +  m2  cos  y)]. 

Using  equations  (1)  and  (2),  the  quantity  in  brack- 
ets in  da  becomes 

h  cos  a  —  h  sin  a  +  cos  y(mi  cos  a  —  m2  sin  a) 
=  cos  be  sin  oq  cos  a  —  cos  50  sin  a  +  cos  y  •  0 
=  cos  b0  sin  (oq  —  a), 

and  in  dS  becomes 

li  sin  a  +  ?2  cos  a  +  cos  y  (mi  sin  a  +  m2  cos  a) 
=  cos  5©  sin  oq  sin  a  +  cos  50  cos  a©  cos  a 

+  cos  y  cos  5 
=  cos  bQ  cos  (a©  —  a)  +  cos  y  cos  5. 

The  corrections  to  be  added  to  the  observed  a  and  8 
of  a  satellite  are  therefore 

/  7  p  (cos  y  —  1 )         «.     •     /  \ 

Ida  =  nr.    . :  cos  5©  sin  (a©  —  a) 

)  2D  sin  y  cos  5 

)  j,.        p(cos  y  —  1)    .        .  ,  . 

\db  =  ^— : : — :  (cos  5©  cos  (a©  —  a) 

2D  sin  y  sin  b 

+  cos  y  cos  8),  (11) 

with  p  =  radius  of  satellite, 

D  =  distance  of  satellite  from  the  camera, 

cos  y  is  obtained  from  equation  (3) ,  and 

aQ,  8©  are  right  ascension  and  declination  of  the 


sun  interpolated  from  the  sun's  ephemeris  for  the 

time  of  observation. 
Because  the  corrections  in  equations  (11)   are  small, 
a  single  entry  in  the  ephemeris  for  the  middle  of  the 
observation  period  will  be  sufficient,  with  a0  and  8Q 
extracted  to  the  nearest  5  seconds  of  time. 

SPECULAR  REFLECTION 

Equations  (11)  apply  to  the  case  of  a  spherical 
satellite  which,  because  of  surface  irregularities,  dif- 
fuses light  in  all  directions  and  hence  from  all  portions 
of  the  illuminated  surface. 

According  to  Snell's  law,  a  parallel  beam  of  light  is 
reflected  in  a  prescribed  direction  from  a  sphere  at 
only  one  point  of  the  surface.  This  point  lies  in  the 
plane  of  the  incident  beam  and  of  the  radius  parallel 
to  the  given  direction  so  that  its  surface  normal,  a 
radius,  also  lies  in  this  plane  and  bisects  the  angle  be- 
tween the  incident  and  reflected  ray. 

In  figure  4,  which  is  similar  to  figure  1,  PF  is  the 
radius  of  the  sphere  which  bisects  the  angle  y  and 
hence  is  also  the  angle  at  F  formed  by  the  incident  ray 
from  the  sun  and  by  the  given  direction  of  the  observer. 
A  distant,  perfectly  reflective,  spherical  satellite  ap- 
pears therefore  as  a  point  source  of  light,  transversely 
displaced  from  its  center  by  a  distance  P"P  — 


P"P  = 


/o      ( i  ~  cos  y\ 

p  sin  y/2  =  Py J  . 


The  corrections  to  a  and  8  are  consequently  the  cor- 
rections in  equations  (11)  of  the  diffusive  case  multi- 
plied by  the  ratio 


"  _  /l  —  cos  y\  ( 1  —  cos  y\  1 
»  ~  \        2        )\        2        / 
_  / 1  —  cos  y \ 
\       2       /    ' 


or  directly 
(da  = 


D  sin  7  cos  5  \        2 

\d*  = 

I  p  (\  -  cos  y\ 

\D  sin  7  sin  8  \        2        / 


(1  —  COS7\  s  .      .      .  x 

I J         cos  50  sin  (oq  —  a) 


(cos  6©  cos  (ocq  —  a) 
+  cos  7  cos  5). 

(12) 


ALTERNATIVE  EQUATIONS  NEAR 
LIMITING  VALUES 

For  y=0°  and  y=180°,  equations  (11)  and  (12) 
are  undefined  as  a  consequence  of  the  implicit  assump- 
tion in  the  development  of  equation  (4)  that  1  and  m, 
the  directions  to  the  sun  and  to  the  earth,  are  independ- 
ent (nonparallel)  vectors.  These  limiting  values  of  the 
phase  are,  however,  never  reached  in  passive  satellite 
photography  as  demonstrated  below  in  connection  with 
figure  6.  The  range  of  values  of  y  near  0°  and  180°, 
excluded  by  those  considerations,  is  sufficiently  large 
to  eliminate  any  computational  difficulties  that  might 
be  anticipated  in  the  evaluation  of  the  corrections  of 
equations  (11)   and  (12)   near  those  singularities. 

On  the  other  hand,  the  zeros  of  sin  8,  cos  8  in  the 
denominator  of  equations  (11)  and  (12)  are  possible 
singularities  and  require  a  modification  of  equations 
(11)  and  (12)  for  the  corresponding  limiting  values 
of  8.  For  8=90°,  da  in  equations  (11)  and  (12)  has 
a  zero  in  the  denominator  corresponding  to  the  fact 
that  right  ascension  is  meaningless  at  the  pole,  hence 
also  its  increment.  This  leaves  as  the  only  real  difficulty 
the  case  of  8=0  for  which  the  given  formulas  for  dS 
become  meaningless.  The  reason  for  this  breakdown 
is  the  arbitrary  choice  of  omitting  dm3  in  setting  up 
equations  (9),  and  this  is  the  only  component  of  the 
vector  dm  that  affects  dS  when  8=0,  as  is  apparent 
from  equations  (8).  For  8=0  therefore,  as  well  as  for 
values  near  0,  it  will  be  necessary  to  determine  dS 
from  the  third  of  equations  (8), 


dS  = 


dm  3 


p(cos  y—  1) 


cos  5       2D  sin  y  cos  5 


(£3  +  rn3  cos  7) 


or 


Figure  4. — Specular  reflection  from  a  spherical  satellite. 


d8  =  ^^ ^(sin6e  +  cos7sin0).        (11') 

2D  sin  7  cos  5 


This  equation  is  an  optional  alternative  for  the  sec- 
ond of  equations  (11),  computable  for  8=0  and  pref- 
erable in  that  region  of  8.  Except  at  the  limits,  both 
yield  identical  results. 

The  corresponding  equation  for  the  case  of  specular 
reflection,  equations  (12),  is 


dd  = 


D  sin  7  cos  5 


{^f^f 


(sin  5©  +  cos  7  sin  5). 

(12') 

Satellite  triangulation  depends  on  the  astronomer's 
star  catalog  for  basic  data  but,  unlike  metric  astron- 
omy, it  operates  and  computes  in  three-dimensional 
space  with  a  Euclidean  metric,  that  is,  with  a  Cartesian 
coordinate  system.  The  singularities  that  require  spe- 
cial consideration  are,  for  the  most  part,  singularities 
of  the  astronomer's  two-dimensional  curvilinear  coordi- 
nate system  which  would  be  avoided  by  dropping  the 
a,  8  concept  and  spherical  trigonometry  at  the  earliest 
possible  opportunity  and  by  adopting  the  methods  of 
analytic  geometry,  preferably  vector  analysis  and 
matrix  calculus.  Experience  at  the  National  Ocean 
Survey  has  shown  that  a  great  many  computational 
programs  can  be  simplified  by  such  a  break  with  tra- 
ditional concepts,  in  addition  to  giving  a  clearer  pic- 
ture of  the  basically  simple  geometric  concepts  in- 
volved. In  the  case  presented  here,  for  example,  once 
an  expression  for  the  camera  direction  has  been  derived 
in  the  form  of  vector  equations  (2),  there  is  really  no 
need  for  explicit  increments  to  a  and  8.  The  increments 
dmi,  dm,2,  dm3  of  equations  (10),  when  added  to  m  of 
equations  (2),  produce  m-j-rfm  which  is  equal  to 
m(a-|-rf«,    8-|-c?8)  ;    the   trigonometric    functions    in- 


volved are  simpler  than  those  required  for  da  and  dS 
and  are  of  general  application  for  all  values  of  a  and  8. 

Again,  the  zero  in  the  denominator  for  y==0°  and 
180°  makes  the  correction  incomputable  for  these  val- 
ues of  y.  Examination  of  figure  1  shows  that  the  cor- 
responding phases  are  "full  moon"  and  "new  moon," 
respectively.  In  the  first  case,  the  phase  correction 
would  be  zero;  and  in  the  second,  no  photograph  is 
possible. 

The  following  geometric  considerations  establish 
limits  for  the  neighborhoods  of  these  points  (y=0°, 
180°)  within  which  the  passive  satellite  cannot  be 
photographed.  In  figure  5,  the  location  of  the  camera 
is  at  E  and  the  sun  is  on  the  horizon  at  elevation  0, 
as  indicated.  For  a  given  fixed  elevation,  the  satellite 
may  occupy  any  of  the  positions  of  the  circle  S,  Sm,  Sm- 

As  S  assumes  various  positions,  the  angle  y  at  5, 
formed  by  the  directions  from  5  to  £  and  from  S  to 
the  sun,  varies  continuously,  increasing  from  a  min- 
imum equal  to  the  elevation  of  the  satellite  eM(  at  the 
position  Sm  to  a  maximum  at  the  position  Sm  where  y 
is  the  supplement  of  e^.  The  points  Sm,  Sm  are  in  the 
vertical  plane  of  the  satellite  containing  the  direction 
to  the  sun,  that  is,  when  the  azimuths  of  the  sun  o©  and 
of  the  satellite  a  differ  by  0°  and  180°,  respectively.  If 
the  angle  of  elevation  or  depression  e©  of  the  sun  is 
different  from  zero,  but  still  fixed,  the  same  argument 
applies  except  that  for  every  position  of  the  satellite 
the  corresponding  angle  is  y+e©.  We  conclude,  there- 
fore, that  to  determine  the  extremal  values  of  y  for 
given  elevations  of  the  sun  and  satellite,  it  is  sufficient, 
as  well  as  necessary,  to  consider  the  situation  when 
sun  and  satellite  are  in  the  same  vertical  plane— the 
plane  of  figure  6. 


0  -* 


[The  length  of  the  generator  SE 
is  unity.     In  right    ASPE: 
PS  =  cos^L  PSE  =cosT.        cosThas 
maximum  when  S  is  at  Sm.J 


Figure  5. — Cone  of  constant  satellite  elevation. 
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Figure  6. — Minimal  phase  angle  for  passive  satellite. 


In  figure  6,  the  camera  is  again  at  E  on  the  surface 
of  the  terrestrial  sphere  with  unit  radius.  The  sun  has 
an  assumed  minimal  angle  of  depression  e©  to  meet 
the  physical  conditions  of  the  problem  and  the  satellite 
moves  in  a  circular  orbit,  concentric  with  the  earth  and 
radius  &>1.  The  sun's  rays  graze  the  earth  at  P,  and 
the  extension  of  this  direction  through  the  point  5  on 
the  satellite  orbit  delineates  the  umbra  for  near-earth 
satellites.  For  such  a  satellite  with  an  elevation  e=ee, 
the  angle  y  would  indeed  be  zero,  but  the  satellite 


would  be  in  the  earth's  shadow  and  hence  excluded 
from  consideration  in  passive  satellite  photography. 
The  satellite  emerges  from  the  shadow  at  5.  A  perpen- 
dicular from  S  onto  a  line  parallel  to  PS  through  E 
creates  the  right  triangle  ERS  so  that  Z  RES  = 
Z  ESP,  which  is  the  minimal  angle  y  for  this  critical 
point  5  of  the  orbit.  The  angle  QOE  =  e0,  and  there- 
fore QE  =  sin  e^,  OQ  =  cos  eQ.  From  right  triangle 
OPS,  we  have  PS  =  ^/¥^T. 


20°r- 


1.2 


1.3 


1.4  1.5  1.6  1.7 

k  (orbit   radius/earth    radius) 

Figure  7. — Minimum  phase  angle. 


1.8 


1.9 


2.0 


Hence,  this  equation.  The  corresponding  elevation  of  the  satel- 

tan  7  =  tan  ARES  lite  at  Point  5™  (of  fig-  5)  wil1  be  e**t  =  eQ+y  (from 

RS             QP             OP  —  00  nS-  6).  The  satellite  remains  visible  and  in  sunlight 

ER  =  qji  _  njjj  =  p<$  _  Q]jj  until  it  reaches  the  horizon  point  H.  At  H,  y  reaches 

or  a  maximum  of  180°  —  e0  which  is  therefore  an  upper 

1  bound  for  the  phase  angle  corrections. 

tan  7  =  — —^^2 •                0-3)  The  graph  of  figure  7  shows  minimum  values  of  y 

vfc2  —  1  —  sine©  derived  from  equation   (13),  with  the  assumed  min- 

The  minimum  value  of  y  for  a  given  angle  of  de-  imum  angles  of  depression  e©  of  the  sun  equaling  10° 

pression  of  the  sun  e0  and  for  a  given  ratio  k  of  the  and  18°,  respectively,  and  the  values  of  k  varying  from 

orbit  radius  to  the  earth's  radius  can  be  computed  from  1.1  to  2. 
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